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An interpretation of the probability flux is given, based on a derivation of its eigenstates and 
relating them to coherent state projections on a quantum wavefunction. An extended definition of 
the flux operator is obtained using coherent states. We present a "processed Husimi" representation, 
which makes decisions using many Husimi projections at each location. The processed Husimi rep- 
resentation reverse engineers or deconstructs the wavefunction, yielding the underlying classical ray 
structure. Our approach makes possible interpreting the dynamics of systems where the probability 
flux is uniformly zero or strongly misleading. The new technique is demonstrated by the calculation 
of particle flow maps of the classical dynamics underlying a quantum wavefunction. 



The probability flux, or probability current, is intro- 
duced in quantum mechanics textbooks as a determinis- 
tic operator that can be calculated, but its connection to 
experiment is often left to the reader's imagination. The 
flux operator, whose expectation over the wavefunction 
gives the traditional flux j (r, p), is denned as 

jr = ^(|r)(r|p + p|r)(r|). (1) 

Here r and p indicate the position and momentum of 
a particle while m is the mass. The concept of "flux 
at a point" seems paradoxical because we say something 
about momentum while also knowing position precisely. 
This raises the question: Can the flux even be measured? 

On the other hand, probability flux vanishes on sta- 
tionary states for systems with time-reversal symmetry. 
This is a shame, since strong semiclassical connections 
between trajectory flow and quantum eigenstates lie com- 
pletely hidden in the universal value of for the flux. 
This letter addresses the problem of the underlying dy- 
namics encoded in stationary states by extending the def- 
inition of the flux to coherent state projections and using 
a "processed Husimi" representation which makes deci- 
sions using many Husimi projections at each location and 
spits out classical rays. Not only does our approach re- 
solve the measurement problem for the flux, we also use 
it as a technique for extracting semiclassical paths from 
a quantum wavefunction even when the flux is zero. 

Several discussions connecting the flux to experimen- 
tal measurement exist in the literature [IHS]; we begin 
our argument instead by identifying the eigenstates of 
the flux operator [4-6J, and present a physical interpreta- 
tion. The flux eigenstates had previously been studied in 
the context of thermal rates [4j [6] and of wavepackets[5]. 
By providing a clear derivation and a novel application 
of the flux eigenstates, we provide a perspective of how 
these techniques can fit into the broader context of wave- 
function analysis. 

We begin by replacing the Dirac basis implicit in Eq. [I] 
with the Gaussian basis defined as 



(v\v ,a) = N^ 2 e-^ 2 ^ 2 , (2) 

where d is the number of dimensions in the system and 

N a = (ay/^) is a normalization constant. There is 
certainly no loss of generality here; the a — >• Dirac delta 
basis limit is always just a step away. Instead, there is 
new capability introduced, as we show below. The Gaus- 
sian basis over all centers is overcomplete, so nothing is 
left hidden from the analysis we develop now. An alter- 
nate derivation using only coherent states confirms the 
validity of our results. Below we are led to a harmonic 
oscillator basis at each location. In the Gaussian basis, 
the modified flux operator is 

j ro ,a = ^ Oo, cr) (r , cr\ p + p |r , a) (r , a\) . 

The eigenstates, projected onto each orthogonal spatial 
dimension z, are obtained using the eigenvalue equation 

jr Q ,a,i \^a,i) = ^a,i \^a,i) •> (3) 

which has a solution of the form 

= |r ,cr) +api |r ,a) . 
Using the following two equations 

(r |p| r , a) = iha' 2 (r - r ) 
and (r , a |p| r ,cr) = 0, we can write 

3r ,*,i \\c,i) = ^ { a (Pi) a l r o,cr) + Pi |ro,cr» • (4) 

Finding the conditions on X a ^ that allow Eq. [4] to be 
written in the form of Eq. [3j we obtain 
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Since (pf ) = we ^ nc ^ ^he vame °f a 
which we obtain the two eigenvalues 
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The eigenstates take the form 

i 

(r| A a>ij ±) = (r| r , a) ± -e* • (r - r ) (r| r , a) , (6) 

(7 

where e$ is the unit vector along spatial direction i. Eq.[6] 
is a linear combination of two functions: the Gaussian 
(Eq. |2| and its derivative. Projection of a wavefunction 
onto the first term can be interpreted as measuring its 
probability amplitude at point r , and projection onto 
second term as measuring its derivative along the i th di- 
rection at the point tq. 

To determine the expectation value of the flux opera- 
tor, we begin by labeling the excited states of the har- 
monic oscillator at position r oriented along the i th di- 
rection 



(r|0) 
(r|l) 

(r|2) 



(r|r ,cr) 
e» • (r - r ) 
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1 (r|r ,cr) ; etc. 



These states form a complete set in which the flux oper- 
ator can be explicitly expressed as the Hermitian matrix 
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where X = X„ 



There are additional sets of 



V,i,+ 4mcr - 

harmonic oscillator states centered at points other than 
ro also with zero components in the flux matrix. 

The complete set of eigenstates |Ai) , | A2) , | A3) , . . . of 
the flux operator expressed in terms of excited states of 
the harmonic oscillator are 



|Ai),|A 2 ),|A 3 ). 



with eigenvalues —A, A, and 0. Measurement by the flux 
operator collapses the wavefunction onto one of these 
eigenstates, the infinite majority of which are in the de- 
generate zero-eigenvalue subspace spanning all excited 
states of the harmonic oscillator above |1). Only the first 
two eigenstates yield non-zero flux values. In the a —> 
limit, the eigenvalues of these two states tend towards 
positive and negative infinity. 
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When expanding the flux expectation value, we can 
use the complete eigenbasis to show that 



Jr ,a,i 



1>) = (V> 



Jr ,o 




A|(V|A 1 )| 2 -A|(^A 2 )r 



(7) 



From Eq. |6j it can be shown that some contributions 
from |(^| 0)| 2 and 1(^1 1)| 2 cancel themselves due to the 
opposite sign of the eigenvalues, and only the cross-term 
(^|0)*(^|1) - (^|0)(^|1)* remains. This form is di- 
rectly related to the following expression for the flux at 
point r 



jr (*(r)) 



2mi 



(^*(r )V^(r ) - ^(r )V^*(r )) . 



Eq. [7| can be rewritten as 

A ih 



= 4^2^ I e i * ( r - r o) |r ,cr) (^|r ,cr) 
-(^\ei • (r-r ) |r ,a)* (^|r ,cr)]. (8) 



The traditional flux operator corresponds to the limit 
a — >• 0, at which point the two terms in Eq.|6]become the 
delta function and its derivative, while the flux values of 
the first two eigenstates become 



lim X a4i ± 



ibOO. 



In addition, there are an infinite number of other eigen- 
states with zero eigenvalues. The infinite value for the 
two non-zero eigenstates is corraborated by the observa- 
tion of Park and Light [4J, however, in their work they 
reach this conclusion in the limit of an infinitely large 
basis |6J. 

A single application of the flux operator at a particular 
point in space almost always results in zero, but occasion- 
ally in an immensely positive or negative value. It is thus 
necessary to perform the averaging over an infinite num- 
ber of measurements to obtain an expression equivalent 
to the textbook flux. 

The prefactors before the Gaussian states in Eq. [6] are 
related to the Taylor expansion 



3 ±^e r (r-r ) 



1 ± -e* • (r - r ). 



This suggests a connection between the flux eigenstates 
and the coherent state, defined as 



(r| r , k , (J) = N d/2 e -{r-r f/^^ .r^ 



(9) 



which is a Gaussian envelope over the plane wave e* k °' r . 
Observing that the phase e zk °' r ° is arbitrary, we Taylor 
expand in the limit of ko<r <C 1 to produce 

(r|r ,k ,cr) w N d/2 e -{r-r ? (1 + . (r _ Tq)) 
w ( r| r , a) + ik ■ (r - r ) ( r| r , a) .(10) 
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Note that the dispersion relation for the free-particle con- 
tinuum is a circle with radius fco = ^ 2 ™ E , which does 
not depend on the orientation of ko. The second terms 
in Eq. 10 and Eq. [6] are proportional to each other when 
ko points along the i th direction. This similarity allows 
us to relate the flux expectation value from Eqs. [7] and [8] 
to coherent state projections as 



lim ( ip 



hko 
Amcj 2 



[|(^|r , k^e h a)Y 



|(^|r ,-/c e;,cr)| 2 ], (11) 



where the traditional flux vector is constructed from the 
components in each direction. 

The representation of quantum mechanical systems 
in terms of phase-space distribution functions, such as 
Wigner's quasi-probability distribution function 018] and 
the Husimi distribution function [9j, allows expressing 
quantum mechanical results into a form that resembles 
classical mechanics [T0HT2] . The standard Husimi projec- 
tion technique consists of an idealized measurement of 
position and momentum consistent with the uncertainty 
principle. The Husimi distribution consists of the expec- 
tation value of the projector onto a coherent state and is 
a phase space probability density. 

Based on the coherent state projection of the flux op- 
erator, we can now develop a deeper intuition. The pro- 
cessed Husimi technique uses coherent states to produce 
maps of the flux that display information about the local 
phase space distribution of the wavefunction beginning 
in real space. As outlined below, the current flow maps 
adhere to the uncertainty principle constraints. Applying 
the processed Husimi representation to a quantum wave- 
function can reverse engineer or deconstruct it, yielding 
the underlying classical rays even when those rays cross 
in several directions, as we demonstrate below. 

Using the uncertainty relation Ax ex 1/Afc ex a and 
setting a — » results in coherent state measurements 
with infinite uncertainty in /c-space and zero uncertainty 
in real space. The traditional flux therefore operates in 
the limit of infinite momentum uncertainty. For small a, 
due to the large momentum uncertainty, coherent state 
projections merely reproduce the probability amplitude 
\i/j(r) | 2 in all directions of k . The flux emerges as a small 
residual which can be retrieved by summing each coherent 
state projection weighted by ko, such that 
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Note that in the limit a — >• 0, the contributing points 
in the /c-space integral reduce to just the orthogonal di- 
rections. The absence of flux in time-reversal symmetric 
states can be interpreted as the mutual cancellation of co- 
herent state projections along each direction in /c-space. 

For larger a, reduced momentum uncertainty allows for 
substantial variation in the coherent state projections be- 
tween different directions of ko. In this regime, we can 
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FIG. 1: States of the circular billiard exhibit defined angu- 
lar momentum and rotational symmetry. The flux (a) is uni- 
formly zero due to time-reversal symmetry. The wavefunction 

(b) , processed Husimi flow (c), and corresponding classical 
paths (d) are shown for an eigenstate of the circular well in 
which the angular and radial components of the wavefunction 
are nearly equal. The classical path in (d) is reproduced in 

(c) to highlight its correspondence with the processed Husimi 
flow. 



use coherent states (Husimi projections) to produce a 
map of the local phase space of a wavefunction. By tak- 
ing snapshots of the phase space at many points across 
a system for larger a, we can process the result to pro- 
duce a semiclassical map showing the dominant classical 
paths contributing to a given wavefunction. Thus, the 
term processed Husimi for these visualizations. Like the 
traditional flux map, processed Husimi flows can be inte- 
grated over lines and surfaces to reveal net current flow. 

Although the calculations reported here are for very 
simple models, we believe that the technique is generic 
and can be applied to virtually any quantum wavefunc- 
tion and model amenable to semiclassical analysis [T3HT5] . 
We begin by demonstrating the processed Husimi tech- 
nique on a circular billiard, which, due to time-reversal 
symmetry, has zero flux. The Schrodinger equation for 
this system can be written in radial form as 



d 2 R(r) 1 dR(r) 
dr 2 r dr 



m 



R(r) = 0. 



Solutions to this equation are simultaneous eigenstates 
of energy and angular momentum, and thus possess the 
good quantum numbers n (number of nodes in the radial 
direction) and m (number of angular nodes). Fig. [TJd 
shows one such state with n w m, which corresponds to 
classical paths that bounce off the boundary at a consis- 
tent 30° angle (see schematic in Fig. [TJi) . 

Because the flux is uniformly zero, it is no help in un- 
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FIG. 2: States of the circular billiard with an applied mag- 
netic field correspond to circular classical paths with a unique 
cyclotron radius. The flux (a), wavefunction (b), processed 
Husimi flow (c), and corresponding classical paths (d) are 
shown for an eigenstate of the circular well with strong mag- 
netic field perpendicular to the plane. The cyclotron radius 
for this state is approximately one-third of the system radius, 
which is corroborated in the processed Husimi flow, but not 
the flux. 



derstanding the dynamics of this system. In Fig. [T]: we 
present the processed Husimi flow of this wavefunction by 
sampling coherent state projections on a regularly-spaced 
grid across the system. With the solid line, we show that 
these projections align perfectly with one of the classical 
paths indicated in Fig. [T]i. In addition, each point in the 
Husimi map contains an additional set of Husimi vectors, 
which do not align with the path. Given that any state 
of a circularly symmetric system must correspond with 
infinitely many classical trajectories related by rotation, 
the processed Husimi at a particular point must reflect 
all rotated paths that intersect there. We indicate one 
classical path and its rotated version in grey in Fig. [T]i. 
The "cross-hatching" pattern arises because two rotated 
classical paths intersect at any point, explaining the sim- 
ilar cross-hatching nodal patterns in the wavefunction. 

Magnetic systems can also be analyzed using processed 
Husimi projections, thus shining light on recent work ex- 
amining flux vortices in quantum dots[l6j H7] , Time- 
reversal symmetry in the circular billiard is broken by the 
magnetic field. To properly represent these states, both 
the momentum operator in the flux operator (Eq. fll) and 
the momentum term iko • 1*0 in the coherent state (Eq. [9| 
must be modified to reflect the canonical transformation 
p — » p — gA/c, where A is the magnetic potential. 

For strong enough magnetic fields, classical trajectories 
are circular with radii corresponding to the cyclotron ra- 
dius. However, the dynamics due to the cyclotron radius 
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FIG. 3: Processed Husimi reconstruction of the classical tra- 
jectories. Coherent state projection vectors for 32 equally- 
space points in /c-space are shown in grey (a-b) for the double 
plane waves (c) defined in Eq. [l3] Because of momentum un- 
certainty, there is spread in the vectors. The flux operator 
(a, blue) averages over the vectors of the coherent state pro- 
jections, while the processed Husimi (b, red) recovers both 
underlying directions. The wavefunction in (c) is representa- 
tive of the areas circled in Fig. [2] 

are not obvious from the wavefunction in Fig. (2]d. Worse 
still, the flux map in Fig. [2^1 (obtained by sampling at 
regularly-spaced points in a grid similar to the one used 
for the processed Husimi flow) appears to indicate two 
circular paths but with strongly different radii. 

The processed Husimi flow in Fig. dramatically clar- 
ifies the dynamics by indicating classical paths consis- 
tent with a cyclotron radius approximately one-third of 
the radius of the system. The cyclotron orbits are offset 
from the center of the system; due to circular symmetry, 
there exist infinitely many such cyclotron orbits rotated 
around the system center. We show a subset of these 
orbits in the schematic in Fig. [2]i, which closely paral- 
lels the processed Husimi flow and exhibits the correct 
cyclotron radius. 

Why does the flux map fail to show the dominant clas- 
sical paths? In Fig. [3j we provide magnified views from 
the full set of coherent state projections, the flux map, 
and the processed Husimi flow corresponding to the cir- 
cles in Figs. [2j We can model this point in the wavefunc- 
tion according to the pure momentum state 

(r) = e ikl,r + e ik2 - r , (13) 

where ki and k 2 are indicated by the arrows in Fig. [3]:. 
Because of the uncertainty of each coherent state projec- 
tion, the full set of projections exhibit a finite spread 
around the generating wavevectors. The processed 
Husimi technique (red) retrieves two independent trajec- 
tories at this point, while the flux (blue) averages over 
them. The flux map amounts to summing the Husimi 
vectors, giving a total flow at each point that does not 
always correspond to the semiclassical dynamics of the 
system. With the processed Husimi flow, we now have 
a complementary representation for revealing the flow 
structure present in the wavefunction. 

Processed Husimi maps also have implications for ex- 
periments measured in a fashion similar to angle-resolved 
photoemission spectroscopy (for a review, see [18J). In 
the ARPES setup, a focused photon beam on a sample 
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emits electrons from the valence band. The energy of the 
photo-emitted electrons incorporates both their bonding 
energies, which can be averaged over, and their kinetic 
energy, which depends on the angle of the beam with 
respect to the sample surface. 

The ARPES response function behaves similarly to co- 
herent state projections with ko proportional to the beam 
angle. By rotating the beam angle around the same point 
of intersection, the response in different directions pro- 
vides the momentum distribution of the wavefunction at 
that point. Perturbations from the known dispersion re- 
lation can then be inserted into Eq. [12] to obtain the flux 
expectation value. 

While a narrow beam could measure the flux vector 
at the intersection point, it will be difficult to distin- 
guish the occasional large perturbation measurements 
from noise. However, wider beams could capture addi- 
tional terms from the Taylor expansion of the coherent 
state in Eq. [lOj producing more reliable measurements. 



Applying the technique at many points across the sample 
would then provide the processed Husimi map. 

We have provided a new interpretation of the flux op- 
erator from the perspective of its eigenstates while con- 
necting them to coherent state projections at the limit 
of infinitesimal spatial spread. Away from this limit, we 
can use coherent state projections to provide a map of 
the classical dynamics underlying a quantum wavefunc- 
tion, permitting us to describe flow even in stationary 
states with zero flux. In systems with magnetic fields, 
we have shown that the flux maps actually correspond to 
the aggregate of such classical flows, which we are able 
to retrieve from the processed Husimi projections. 
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